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C ■ ABSTRACT 

1—5 ; 

The efficiency of thin disk accretion onto black holes depends on the inner 
boundary condition, specifically the torque applied to the disk at the last stable 
Kj^ , orbit. This is usually assumed to vanish. I estimate the torque on a magnetized 

cn [ disk using a steady magnetohydrodynamic inflow model originally developed by 

^ ! [I'akahashi et al. 1990| . I find that the efficiency e can depart significantly from 



c^ 



O, 

Q\ , the black hole. 

o 



the classical thin disk value. In some cases e > 1, i.e. energy is extracted from 
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1. Introduction 



>< 
V^ ■ The dynamics of accretion close the event horizon of black holes is of considerable 

interest to astronomers because that is where most of the accretion energy is released and 

because strong-field gravitational effects may be evident there. 

Much of the work on thin disks around black holes assumes that the disk ends at or near 
r = Tmso, the radius of the marginally stable circular orbit. Inside Vmso, material is assumed 
to nearly free-fall onto the black hole (e.g. [Bardeen 1970| , [I'horne 1974| , [Abramowicz et aT 



1978|) . This implies the so-called "no-torque" boundary condition on the disk, and leads 



to a disk whose surface brightness vanishes at the inner edge QShakura fc Sunyaev 1973|) . 



The no-torque boundary condition leads to the classical estimate for the accretion efficiency 
eo = 1 — £{rmso)/c^, where S{rmso) = —Utirmso) is the "energy at infinity" per unit mass 



^also Harvard-Smithsonian Center for Astrophysics, MS-42, 60 Garden St., Cambridge, MA 02138, USA 
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of a particle at the marginally stable orbit, eo ranges between 1 — 2^/2/3 ~ 0.057 for a 
nonrotating (a = 0) black hole, to 1 — 1/a/3 ~ 0.42 for a prograde disk around a maximally 
rotating (a = 1) black hole. 

It was realized early on, however, that magnetic fields might alter the dynamics of 
the accreting material and hence the accretion efficiency (see, e.g., J. M. Bardeen quoted 
in fThorne 1974| ). More recently, it has been argued that the magnetic fields threading the 



accreting material should be strong enough to be dynamically important, but not so strong 
that the fields can be regarded as force- free ( [Krolik 1999| ). 

In this Letter I estimate the magnetic torque on the inner edge of a thin accretion disk 
in the Kerr metric. I use a model for the infiow that is based on earlier work by p?akahash] 



et al. 19901 , following Phinney 1983| and Camenzind 1986a,1986b,1989. 



2. Model 

Consider a thin disk in the equatorial plane of the Kerr metric. Because the disk is 
thin, c^s/(? <C 1, so the relativistic enthalpy r/ ^ 1. The disk is magnetically turbulent 
( Palbus fc Hawley 1991| ) with magnetic energy density B'^/{8tt) ^ 2apcl (|Hawley, Gammie 



I&: Balbus 1995|) , where a is the usual dimensionless viscosity of accretion disk theory. 
I will assume that a <C 1 so that magnetic fields make a negligible contribution to the 
hydrostatic equilibrium of the disk. Then the disk has an inner surface at r = rj„, with 
rmso — fin ~ Tmsoics/cY^ SO rj„ ^ Tmso- At the inner surface of the disk I imagine that 
magnetic field lines rise approximately radially through the disk atmosphere and force the 
atmosphere to corotate with the disk surface. High enough in this corotating atmosphere, 
the gradient of the effective potential changes sign and the gas begins to stream inward 
along field lines toward the black hole. As gas leaves the disk, the character of the fiow 
changes, becoming less turbulent and more nearly laminar. 

This picture leads one to consider a steady, axisymmetric infiow close to the equatorial 
plane of the Kerr metric. In what follows, I will work exclusively in Boyer-Lindquist 
coordinate t, r, 6', 0, and follow the notational conventions of [Misner, I'horne, fc Wheeler| 



1973i I assume that, in the poloidal plane, the fiuid velocity and magnetic field are 
purely radial, i.e. that m^ = and F^^ = 0, where u^ is the four-velocity and F^^, is the 
electromagnetic field ("Maxwell") tensor; recall that in the nonrelativistic limit, denoted 
— i>, Fr(j){6 = vr/2) -^ rBg, where B is the magnetic field three-vector. I also assume that all 
fiow variables are functions only of r, i.e., that de = 0. This one- dimensional fiow is similar 
in spirit to the old [Weber fc Davis 1967| model for the solar wind, turned inside-out. As in 
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the Weber-Davis model, the magnetic field has a monopolar geometry. 

I will also assume perfect conductivity, so that the electric field vanishes in a frame 
comoving with the fluid: 

u^F^, = 0. (1) 

Together with the symmetry conditions, this leaves F^j, with six nonzero components: 
Fte {-^ —f^Eg), Fro (^ —rB^), and Fg^ (-^ —r'^Bj.), and gives the relation 
Fre = —{Fteu^ — F0^u'^)/u^ . Thus the electromagnetic field has only two degrees of 
freedom. 

The flow symmetries reduce Maxwell's equations to 

drFte = (2) 

and 

drFe^ = 0. (3) 

The flrst equation is the relativistic "isorotation" law, which can be rewritten Fte = ^pFg^, 
where Qp is the rotation frequency u'^/u^ at the radius where u^ = 0. The second equation 
is the relativistic equivalent of V ■ -B = 0. 

Conservation of particle number leads to a conserved "rest mass flux" , per unit 6, 

Fm = 2'Kr'^pu'' -^ 27ir'^pvr, (4) 

here p = particle number density multiplied by the rest mass per particle and v is the 
three- velocity. The conserved angular momentum flux is 

F. = 2.rX' = 2.r^ ^«' - ^g^) (5) 

3 B4,Br, 



2-Kr [pVrVfj, 



An 



where T^^^y is the stress-energy tensor, V = 1 — 2r^/r + a^r'^/r'^, Vg = GM/c^. The conserved 
mass-energy flux is 

-2nr'T[ = 2nr' {-putu^ - ^^^) (6) 
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and, consistent with the thin disk approximation, I have neglected the thermal energy's 
contribution to Fe- Notice that the accretion efficiency is given by e = 1 — Fe/ Fm- 
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The final equation needed to close the system is the normalization of four-velocity, 

w% = -c\ (7) 

Henceforth I will set Vg = c = —Fm = 1-0 

Equations (|l]) through (^ describe a one dimensional inflow model with a series of 
algebraic relations. This model is physically identical to that developed by P?akahashi et al 



|1990| , although cast is somewhat different notation. Also, unlike [Takahashi et al. 1990| , I 



specialize to the case where the inflow is anchored in a thin disk at the marginally stable 
orbit. This sets the inflow in a specific astrophysical context and allows one to relate it 
to an accretion efficiency. These inflow solutions have also been explored by pamenzi"nd 



1994| , [Uamenzind 1996| and, in unpublished work, by Camenzind and Englmaier. 



3. Boundary Conditions 

The inflow model has 6 dynamical variables: p, m*, u''', u'^, Fto, and F^^, and 6 conserved 
quantities: Fm, F^, Fe, Fg^f^Qp, and u^m^, so the equations are fully integrated. What sets 
the conserved quantities? 

The mass flux Fm has been normalized to —1, but more physically it is set by conditions 
in the disk at large radius. The normalization of four- velocity gives m'^m^ = — 1. F^^ is 
related to the magnetic flux emerging from the inner edge of the disk. This is presumably 
determined by the action of magnetohydrodynamic (MHD) turbulence in the disk and the 
interaction of the disk with the inflow itself. I will treat it as a parameter of the problem, 
and estimate reasonable values for it later. Three quantities remain to be specified. 

It seems reasonable to require that the energy and angular momentum fluxes be 
continuous across the boundary between the inflow and the disk. It also seem reasonable 
to require that the inflow four-velocity match continuously onto the disk. Happily, these 
requirements do not overconstrain the problem. If Qp = i7(rj„) and u'^{rin) = 0, then one 
can show that Fp = flFp + {S — IVL)Fm at rj„; here I = u^. This ehminates two degrees 
of freedom by expressing Fp in terms of Fi and by fixing VLp. The energy and angular 
momentum flux are then continuous if the viscous (more properly, turbulent) angular 
momentum flux of the disk just outside r^n matches the electromagnetic angular momentum 



^Physical units may be recovered as follows. Length: r^; time; Tg/c] mass: —FMc/rg\ Fg, 



GM{^Fm)^/^c-^/^, while Fm ~ -Mr/{2H). 
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flux of the inflow just inside rj„: 

T;(visc) = r;(EM). (8) 

This condition will be automatically satisfied once the disk relaxes to a steady state, 
since the disk will increase its surface density until it can carry off the entire outward 
electromagnetic angular momentum (and energy) flux from the inflow. 

The final degree of freedom {F^, or equivalently, Fe) is fixed by the condition that 
the flow pass smoothly through the fast critical point. Thus Fl and Fe emerge as 
"eigenvalues" of the solution. 

4. Solutions 

Once the model parameters are fixed, the resulting set of nonlinear algebraic equations 
must be solved numerically. I obtain Fl, and the location of the fast critical point {rf,u^f) 
via simultaneous solution of 

durFE{rf,u},FL) = 0, (9) 

drFE{rf,u},FL) = 0, (10) 

FE{rf,u},FL)-FEir^so,0,FL) = 0. (11) 

Here I have used the fact that the critical point is a saddle point of ^^(r, m^; Fl). I use 
the multidimensional Newton- Raphson method of [Press et al. 1992| , and for simplicity I 
evaluate the derivatives numerically. 

There is one subtlety involved in the solution. In calculating FE{r,u'^,FL) one must 
solve a quadratic equation for u* (or equivalently m*^), so one must decide which root is 
physical. In general this is a nontrivial matter, since the physical solution skirts a region 
in the r, u^ plane where the discriminant of the quadratic vanishes (the solution makes a 
smooth transition from one branch of the solution to another). Fortunately it turns out that 
one branch is appropriate in the neighborhood of Vmso and the other in the neighborhood of 
the fast point, so in practice this subtlety is easily dealt with. 

Figure 1 shows u'^{r) for a solution with a = 0.95 and Fg^ = 6. The dot marks the 
fast point at (r,^'^) = (1.37,-0.26). This solution has Fe ^ 0.04 and Fl ^ 1.23, so 



^The Alfven point does not impose any new condition on the flow, since all trans- Alfvenic solutions pass 



smoothly through the Alfven point; see, e.g., Phinncy 1983. The slow point is absent because the flow is 
cold. 
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e = 1 + Fe ~ 1.04, i.e. energy is extracted from the black hole. Evidently magnetic fields 
can significantly alter the efficiency from the classical value. 

Figure 2 shows contours of the eigenvalue Fl in a survey over the a — Fg^ plane (see 
also Fig. 3 of |Camenzind 1994] ) . For strong field and large black hole spin (up and to 



the right of the heavy solid line in the figure), F^ > 0. [Takahashi et al. 1990| have shown 



that Fi > if the field rotation frequency Vtp is exceeded by the characteristic rotation 
frequency of the space time, 2a/(r^(l + a? /r"^ + 2a^/r^)), at the Alfven point. 

Figure 3 shows contours of the accretion efficiency e evaluated over a portion of the 
a — Fq^ plane (see also Fig. 4 of [Uamenzind 1994] ). The contours in Figure 3 are located 



at intervals of Ae = 0.05. For strong field and large black hole spin (up and to the right of 
the heavy solid line in the figure), Fe > 0, that is, energy is extracted from the black hole. 
As P?akahashi et al. 1990| have shown, this can only happen if the Alfven point lies within 
the ergosphere (r < 2). The heavy dashed line corresponds to the classical efficiency of 
prograde disk around a maximally rotating black hole, e ~ 0.42. 



5. Astrophysical Discussion 

What is an astrophysically sensible value for the crucial magnetic field strength 
parameter Fg^7 I will make a purely Newtonian estimate, for clarity and because the 
relativistic corrections are likely to be smaller than the other sources of uncertainty. 
Suppose that the radial field leaving the disk is = fB^, where subscript d denotes a quantity 
evaluated in the disk and / ^ 1. It is a result of numerical models of disk turbulence that 
i?J/(87r) ^ 2apdC^ad (^-S- [Hawley, Gammie fc Balbus 1995] ). Using the usual steady state 



disk equation 37rSz/ ^ M ^ -2FM{H/r), I find 

Fe^ ^ r^Br ^ 2.3/r;^/^ (12) 

If the inner edge of the disk is at r^ = 6, F^,^ ^ 8-8/5 so the region of parameter space 
shown in Figures 2 and 3 is likely relevant to disks. 

The infiow model thus suggests that the presence of a modest magnetic field, and 
the accompanying torque on the inner edge of the disk, can significantly increase the 
efficiency of thin disk accretion onto black holes. For a = and Fg^ -C 1, e ~ eo + O.OljFe^l, 
while for Fg^ = 4, e = 0.165. The added luminosity manifests itself as an increase in 
the surface brightness of the disk due to the torque applied at its inner boundary. In the 
same quasi- Newtonian spirit as the estimate of Fg^, the surface brightness profile becomes 
oc 1 — /5(rmso/r)^/^ (see ^hapiro fc 'leukolsky 1983| ) and /3 ;^ 1 in the Shakura-Sunyaev 
model. Here /3 = 1 - Fi(EM)//„,, ^ 1. 
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The inflow model also suggests that, because material accretes with a smaller specific 
angular momentum than it would in the absence of magnetic fields, it is more difficult to 
manufacture a rapidly spinning hole by disk accretion. An equilibrium spin is reached when 
Fl = 2aFE (e.g. jPopham fc Gammie 1998| , eq. 10). For Fq^ = 4 this equilibrium value is 
reached at a ^ 0.7. 



A major limitation of this study (which should be regarded as an instructive example 
rather than a source for estimating efficiencies), is that I have ignored the vertical structure 
of the infiow. Crudely speaking, one might expect that field lines which do not lie in 
the midplane are more lightly loaded {\Fg^\ is larger) so that specific energy and angular 
momentum fiuxes might increase away from the midplane, until at sufficient latitude one 
reaches a field line which would rather infiate away than remain tied between the infiow 
and the disk. At high latitude, then, the outward electromagnetic energy fiux might emerge 
in the form of a wind, and be better described by the force-free magnetosphere model of 
Blandford fc Znajek 1977 . 



Another major limitation is the simplified, monopolar field geometry. This limitation 
can be overcome by direct numerical integration of the basic equations. Evidently numerical 
models of infiow inside the marginally stable orbit would be enormously interesting. 
Fortunately they are now practical, at least within the MHD approximation. 

Stu Shapiro's thoughtful advice has greatly improved this paper. I am also grateful to 
Eric Agol, Roger Blandford, Jeremy Goodman, John Hawley, Julian Krolik, Jochen Peitz, 
and Gordon Ogilvie for their comments, and to Phil Myers and the Radio Group at the 
Center for Astrophysics for their hospitality. 
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Fig. 1. — An example solution with a = 0.95 and Fg^f, = 6. The upper left panel shows u^{r), 
the upper right / = U(j,{r). The dot shows the location of the fast critical point. The lower left 
panel shows the electromagnetic energy density {T^y{EM)u^u^; solid line) and the rest-mass 
energy density of the accreting particles (p; dashed line) measured in a frame comoving with 
the fluid. The lower right panel shows the electromagnetic, particle, and total energy fluxes 
as a function of radius. 
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Fig. 2. — Contours of constant angular momentum flux F^ in the a, Fq^ plane. The zero 
angular momentum flux contour is marked as a heavy solid line; above and to the right of 
this line Fl > 0. The light contours are linearly spaced at AFl = 1. 
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Fig. 3. — Contours of constant efficiency e = 1 + Fe in the a, Fg^ plane. The unit efficiency 
contour is marked as a heavy sohd hne; above and to the right of this hne Fe > 0, i.e. energy 
is being extracted from the black hole. The light contours are linearly spaced at Ae = 0.1. 
The heavy dashed contours lies at e ?« 0.42, the classical efficiency of prograde disk around 
a maximally rotating black hole. 



